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1 Introduction 

Traditionally, condensed matter physics has focused on the investigation of 
perfect crystals. However, real materials usually contain impurities, disloca- 
tions or other defects, which distort the crystal. If the deviations from the 
perfect crystalline structure are large enough, one speaks of disordered sys- 
tems. The Anderson model [1] is widely used to investigate the phenomenon of 
localisation of electronic states in disordered materials and electronic trans- 
port properties in mesoscopic devices in general. Especially the occurrence 
of a quantum phase transition driven by disorder from an insulating phase, 
where all states are localised, to a metallic phase with extended states, has 
lead to extensive analytical and numerical investigations of the critical prop- 
erties of this metal- insulator transition (MIT) [2-4] . The investigation of the 
behaviour close to the MIT is supported by the one-parameter scaling hypoth- 
esis [5,6]. This scaling theory originally formulated for the conductance plays 
a crucial role in understanding the MIT [7]. It is based on an ansatz inter- 
polating between metallic and insulating regimes [8]. So far, scaling has been 
demonstrated to an astonishing degree of accuracy by numerical studies of the 
Anderson model [9-13]. However, most studies focused on scaling of the local- 
isation length and the conductivity at the disorder-driven MIT in the vicinity 
of the band centre [9, 14, 15]. Assuming a power-law form for the d.c. conduc- 
tivity, as it is expected from the one-parameter scaling theory, Villagonzalo et 
al. [6] have used the Chester- Thelhmg-Kubo-Greenwood formalism to calcu- 
late the temperature dependence of the thermoelectric properties numerically 
and showed that all thermoelectric quantities follow single-parameter scaling 
laws [16,17]. 

In this chapter we will investigate whether the scaling assumptions made 
in previous studies for the transition at energies outside the band centre can be 
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reconfirmed in numerical calculations, and in particular whether the conduc- 
tivity <7 follows a power law close to the critical energy E c . For this purpose we 
will use the recursive Green's function method [18, 19] to calculate the four- 
terminal conductance of a disordered system for fixed disorder strength at 
temperature T = 0. Applying the finite-size scaling analysis we will compute 
the critical exponent and determine the mobility edge, i.e. the MIT outside 
the band centre. A complementary investigation into the statistics of the en- 
ergy spectrum and the states close to the MIT can be found in Chapter [20]. 
An analysis of the mathematical properties of the so-called binary-alloy or 
Bernoulli- Anderson model is done in Chapter [21]. 

2 The Anderson Model of Localisation and its 
Metal-Insulator Transition 

The Anderson model [1,2] is widely used to investigate the phenomenon of 
localisation of electronic states in disordered materials. It is based upon a 
tight-binding Hamiltonian in site representation 



where |i) is a localized state at site i and Uj are the hopping parameters, which 
are usually restricted to nearest neighbours. The on-site potentials are ran- 
dom numbers, chosen according to some distribution P(e) [22,23]. In what 
follows we take P(e) to be a box distribution over the interval [— W/2, W/2], 
thus W determines the strength of the disorder in the system. Other distri- 
butions have also been considered [2,3,24]. 

For strong enough disorder, W > W C (E = 0), all states are exponen- 
tially localized and the respective wave functions W(r) are proportional to 
c (-|r- r o|/C) f or l ar ge distances |r — ro|. Thus, \P is confined to a region of some 
finite size, which may be described by the so-called localisation length £. In 
this language extended states are characterised by £ — ► oc. Comparing £ 
with the size L of the system one can distinguish between strong and weak 
localisation, for £ <C L and L < £, respectively 3 . Here we also assume that 
the phase-relaxation length 1$ 3> L. Otherwise, the effective system size is 
determined by t$. 

It turns out that the value of the critical disorder strength W c depends on 
the distribution function P(e) and the dimension d of the system. In absence 
of a magnetic field and for d<2 all states are localized 4 , i.e. W c — [7,8]. 

3 We note that the phrase weak localization in the context of the scaling theory 
is often used with a specific meaning, namely the onset of localization in large 2- 
dimensional samples, where the conductance decreases logarithmically with scale 



4 Strictly speaking, this is only true if Tt belongs to the Gaussian orthogonal 
ensemble [25]. 
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For systems with d = 3 the value of W c additionally depends on the Fermi 
energy E and the curve W C (E) separates localized states (W > W C (E)) from 
extended states (W < W C {E)) in the phase diagram [22,23,26]. If instead 
of E the disorder strength is taken as a parameter, there will be a critical 
energy E C (W) — also called the mobility edge — and states with \E\ < E c are 
extended and those with \E\ > E c localized yielding the same phase boundary 
in the (E, W^-plane. At the mobility edge, states are multifractals [27]. The 
separation of localized and extended states is illustrated in Fig. 1, which shows 
a schematic density of states (DOS) of a three-dimensional (3D) Anderson 
model. Since for T = localized states cannot carry any electric current, the 




Fig. 1. Typical DOS of a 3D Anderson model for fixed W < W c . The states in 
the grey regions are localized, otherwise they are extended. The mobility edges are 
indicated at ±E C . Also indicated is the power-law behaviour of cr(E) (dashed lines) 
close to ±E C according to (2). 



system shows insulating behaviour, i.e. the electric conductivity a vanishes 
for \E\ > E c or W > W c . Otherwise the system is metallic. Therefore, the 
transition at the critical point is called a disorder-driven MIT. 

For the MIT in d = 3 it was found that a is described by a power law at 
the critical point [2], 



a(E) = 



Co 



1 - 
0, 



\E\ < E c 
\E\ > E c 



(2) 



with v being the universal critical exponent of the phase transition and Co a 
constant. The value of v has been computed numerically by various methods 
[2,9-11] and was also derived from experiments [28,29]. The results range from 
1 to 1.6, depending on the distribution P(e) and the computational method [3] 
used. 
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Moreover, Wegner [30] was able to show that for non-interacting electrons 
the d.c. conductivity a obeys a general scaling form close to the MIT, 

a(e,u) ^b 2 - d a{b 1/v e 1 b z u J ) . (3) 

Here e denotes the dimensionless distance from the critical point, u> is an 
external parameter such as the frequency or the temperature, b is a scaling 
parameter and z is the dynamical exponent. For non- interacting electrons 
z = d [31]. Assuming a finite conductivity for uo = 0, one obtains from (3) 

a(s,0) (4) 

where s = |1 - £/£ c |- With d = 3 this gives (2). 

3 Computational Method 

An approach to calculate the d.c. conductivity from the Anderson tight- 
binding Hamiltonian (1) is the recursive Green's function method [18,19,23]. 
It yields a recursion scheme for the d.c. conductivity tensor starting from 
the Kubo-Greenwood formula [32]. Moreover, this method allows to com- 
pute the density of states and the localization length as well as the full 
set of thermoelectric kinetic coefficients [33] . Parallel implementations of the 
method are advantageous [34,35]. The method is therefore a companion to the 
more widely used transfer-matrix method [36, 37] or iterative diagonalisation 
schemes [38,39]. 

3.1 Recursive Green's Function Method 

Let TL = Y^ij denote our hermitian tight-binding Hamiltonian. The 

single particle Green's function Q ± {z) is defined as [40] (z* — TC)Q ± = 1 where 
z = E ± Y) is the complex energy and the sign of the small imaginary part 7 
distinguishes between advanced and retarded Green's functions, G~(E — iO) 
and Q + {E + iO), respectively [40]. Equivalently, Q ± can be represented in the 
basis of the functions 

- H %j )G% = 6 l3 , (5) 

where Gfj is the matrix element (ilG^ \j}- We note that for a hermitian Hamil- 
tonian Gr. = (G+f. 

If TL contains only nearest-neighbour hopping matrix elements, (5) can be 
simplified using a block matrix notation. This is equivalent to considering the 
system as being built up of slices or strips for 3D or 2D, respectively, along 
one lattice direction. In what follows all quantities written in bold capitals are 
matrices acting in the subspace of such a slice or strip. For 2D and 3D these 
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Fig. 2. Scheme of the recursive Green's function method for a 3D system. The 
new Green's function G^ JV+1 ' can be calculated from the old Hamiltonian Hjv (light 
grey), the new slice Hamiltonian Hjv+i (dark grey) and the coupling tjv (solid 
arrows) . 



are matrices of size M x M and M 2 x M 2 , respectively, where M is the lateral 
extension of the system (cf. Fig. 2). The left hand side of (5) is then given as 



/ 







-H M _! (zn-Hu) 



H 



H 



i+l,i 



(z±I-H 



i,»+l 
»+l,*+l 



) H 






i+l,»+2 



V 



G t-l,j 

G i+1,3 



(6) 



where i and j now label the slices or strips. From this expression one can easily 
see that (5) is equivalent to 



(z±I - H 4 ,)G± - H^Gt, 



H 



l G i+l,j - 



(7) 



Using the hermiticity of TL we define the hopping matrix tj = H iii+ i (and 
hence tj = H^-i) connecting the zth and the (i + l)st slice. Now, we consider 
adding an additional slice to a system consisting of N slices. The Hamiltonian 
of this larger system can be written as [19] 



H (N+1) _ >H y+t JV +t] v + H JV+ i lJV+1 (i,j<N). 



(8) 



The first and the last terms describe the uncoupled iV-slice and the additional 
1-slice system. Using t n as an "interaction" the Green's function G^ +1 ' of 
the coupled system can be calculated via Dyson's equation [19,40], 



(9) 
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In particular, we have 



I-H 



r (N+l) _ r (JV) 



r (N+l) 
"i.JV+l 



JV+1, JV+1 



ut r Wi 

N NN N 



r"W v J + r , ( w+1 ) +t r" 1 - 

^JJV l N*-*N+l, JV+1 L JV"JVj 



i-t 



r (JV) p (JV+l) 



G 



(JV+1) 
JV+lj 



G 



JV+1, JV+1 



(i < JV) 



(JV+1) 

JV+1, JV+1 u JV"JVj 



tlvG^V (j<iV). 



(i,j<JV) 



(10a) 
(10b) 
(10c) 
(lOd) 



With (10) the Green's function can be obtained iteratively. Additionally, there 
are two kinds of boundary conditions which must be considered: across each 
slice and at the beginning and the end of the stack. The first kind does not 
present any difficulty and usually hard wall or periodic boundary conditions 
are employed. The second kind of boundary is connected to some subtleties 
with attached leads which will be addressed in Section 6. 



3.2 Density of States and D.C. Conductivity 

The DOS is given in terms of Green's function by [40] 

P(E) = -^I m ^+ = -^Imf>G+ (11) 

and the d.c. conductivity a is 

a = -^A Tr [p Im Q+ p Im Q+] . (12) 

Here, fl denotes the volume of the system and m the electron mass. Using 
for the momentum the relation p — ^\H,x] one can rewrite (12) in position 
representation 

2a( n n ~] 

a = hh^ Tl yT, G t * G 7* « - i| E( G ^ - x " > ( 13 ) 

l i,j i ) 

where Xi is the position of the zth slice. 

Starting from these relations and using the iteration scheme (10) one can 
derive recursion formulae to calculate the properties for the (N + l)-slice sys- 
tem. The results are expressed in terms of the following auxiliary matrices 
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(14a) 
(14b) 

(14c) 
(14d) 
(14c) 



The derivation can be simplified assuming the new slice to be at xat+i = 0. 
This leads, however, to corrections for the matrices Bat and because the 
origin of Xi has to be shifted to the position of the current slice in each iteration 
step. The corrections are 



B'at — Bat + iC^ + 10^ + -tJ ir (Rjv — Rjv)t 
= — i— t]y(Rjv — R] v )tjv 



(15a) 
(15b) 



Here we have used the identity 



lJ2 G m G m = ^(Gnn - G NN ) = i^(Rjv - K) = ~^^n ■ (16) 

The derivation of the recursion relations is given in Refs. [18, 19,23,41], it 
yields the following expressions 



s (^+i) =s W+Tr{R Ar+1 (F Ar +I)}, 

s (N+i) = S W +T r{Rc(B Ar R w+1 ) + C+Rj v+1 C w R w+ i} 



Rat+i 
Bat+i 



z \ — H 



Ar+1,AT+1 



tjV +1 RAf+l 



B 



A r 



— tjyRArtAf 



2C A 7Rj v+1 C Ar 



R 



C^+i — tj v+1 RAr+iC A ^Rj v+1 tAT+l , 



tAr+lRAT+l^ArB-Ar+ltAr+l 



Fat+1 — tj v+1 RAT+l(FAr + I)RAT+ltAT+l • 

The DOS and the d.c. conductivity are then given as 



(17a) 
(17b) 

(17c) 

(17d) 

(17e) 
(17f) 
(17g) 
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( N + 1 )(E) = - s (Af+1) (18) 

a {N+1) (E) = — - s (N+1) (19) 

For a comparison with the scaling arguments, we convert the conductivity 
into the two-terminal conductance as 

M 2 

92 = o— (20) 

with L = N + 1. In distinction to the usual use of the recursive scheme which 
constructs a single sample with L ^> M, we shall have to use many different 
cubic samples with L = M. 

We note that it is also possible to calculate the localisation length £(E) by 
the Green's function method. The value of £(E) is connected to the matrix 
r + 

W)=- n^^ ln i TrG - ( <- (21) 

The recursion relation for ^ N+1 \E) is 
1 1 



^( N + 1 )(E) N + l S i 

(22b) 



(N+1) , (22a) 



Tr . r + (N+1) 



4 Finite-Size Scaling 

For finite systems there can be no singularities induced by a phase transition 
and the divergences at the MIT are always rounded off [42]. Fortunately 
the MIT can still be studied using a technique known as finite-size scaling [2] . 
Here we briefly review the main results taking the dimensionless four-terminal 
conductance 34 of a large cubic sample of size L x L x L as an example. We 
note that similar scaling ideas can also be applied to the reduced localisation 
length £/L. In order to obtain 174 of the disordered region only, we have to 
subtract the contact resistance due to the leads. This gives 

94 92 M 

Here N = N{E) is the number of propagating channels at the Fermi energy 
E which is determined by the quantization of wave numbers in transverse 
direction in the leads [43,44]. 

Near the MIT one expects a one-parameter scaling law for the dimension- 
less conductance [7, 30, 42] 
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gi(L,s,b) = T 



(24) 



where b is the scale factor in the renormalisation group, x is a relevant scaling 
variable and v > is the critical exponent. The parameter e measures the 
distance from the the mobility edge E c as in (2). However, recent advances in 
numerical precision have shown that in addition corrections to scaling due to 
the finite sizes of the sample need to be taken into account so that the general 
scaling form is 



g 4 (L,e,b) = T 



p X (e)6 1/v ,0(e)6- B 



(25) 



where <f> is an irrelevant scaling variable and y > is the corresponding ir- 
relevant scaling exponent. The choice b = L leads to the standard scaling 
form 5 

(26) 



with F being related to T. For E close to E c we may expand F up to order 
riR in its first and up to order n\ in its second argument such that 



34 ( 



i(L,e) = 4> n 'L- n ' y F nl (xL llv ) with 

n'=0 

F nl { X L 1 l") = Y j a n > nX n L n/v . 



(27a) 
(27b) 



n=0 



Additionally \ an d 4> may be expanded in terms of the small parameter e up 
to orders tor and m\, respectively. This procedure gives 



TOR 

E 6 - 



(e) 



mi 

E 



(28) 



From (26) and (27) one can see that a finite system size results in a systematic 
shift of g4(L,e — 0) with L, where the direction of the shift depends on the 
boundary conditions [42]. Consequently, the curves gi{L,e) do not necessarily 
intersect at the critical point e = as one would expect from the scaling law 
(24). Neglecting this effect in high precision data will give rise to wrong values 
for the exponents. 

Using a least-squares fit of the numerical data to (27) and (28) allows us to 
extract the critical parameters v and E c with high accuracy. One also obtains 
the finite-size corrections and can subtract these to show the anticipated scal- 
ing behaviour. This finite-size scaling analysis has been successfully applied to 
numerical calculations of the localisation length and the conductance within 
the Anderson model [3, 14]. 



^The choice of b is connected to the iteration of the renormalisation group [42]. 
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5 MIT at E = 0.5 t for varying disorder 
5.1 Scaling of the Conductance 

We first investigate the standard case of varying disorder at a fixed energy [14] . 
We choose Uj ^ for nearest neigbours i, j only, set tij = t and E — 0.5 1 
which is close to the band centre. We impose hard wall boundary conditions 
in the transverse direction. For each combination of disorder strength W and 
system size L we generate an ensemble of 10000 samples. The systems under 
investigation are cubes of size L x L x L for L = 4, 6, 8, 10, 12 and 14. For each 
sample we calculate the DOS p(E, L) and the dimensionless two-terminal con- 
ductance c/2 using the recursive Green's function method explained in Section 
3. Finally we compute the average DOS (p(E,L)), the average conductance 
((^(-E 1 , L)) and the typical conductance exp(ln gi(E, L)). 

The results for the different conductance averages are shown in Figs. 3 and 
4 together with respective fits to the standard scaling form (26). Shown are the 
best fits that we obtained for various choices of the orders of the expansions 
(27, 28). The expansion orders and the results for the critical exponent and 
the critical disorder are given in Table 1. In Fig. 5 we show the same data as in 
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Fig. 3. Average dimensionless conductance vs disorder strength for E = 0.5 1. Sys- 
tem sizes are given in the legend. Errors of one standard deviation are obtained from 
the ensemble average and are smaller than the symbol sizes. Also shown (solid lines) 
are fits to (26) for L = 8, 10, 12 and 14. 



Figs. 3 and 4 after the corrections to scaling have been subtracted indicating 
that the data points for different system sizes fall onto a common curve with 
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Fig. 4. Logarithm of the typical dimensionless conductance vs disorder strength for 
E = 0.5 1. System sizes are given in the legend. Errors of one standard deviation 
are obtained from the ensemble average and are smaller than the symbol sizes. The 
solid lines are fits to (26) . 



average W min /t VK max /i 


n-R. 
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TOr TOi 
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(g 4 ) 15.0 18.0 


2 
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1.55 ±0.11 16.47 ±0.06 




(ln ff4 ) 15.0 18.0 


3 


1 


1 


1.55 ±0.18 16.8 ±0.3 O.i 


± 1.0 



Table 1. Best-fit estimates of the critical exponent and the critical disorder for 
both averages of g 4 using (26). The system sizes used were L = 8, 10, 12, 14 and 
L — 4,6, 8, 10, 12, 14 for (g) and (In g), respectively. For each combination of disorder 
strength W and system size L we generate an ensemble of 10000 samples. 
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two branches as it is expected from the one-parameter scaling theory. The 
results for the conductance averages and also the critical values are in good 
agreement with transfer-matrix calculations [9, 14]. 

5.2 Disorder Dependence of the Density of States 

The Green's function method enables us to compute the DOS of the disordered 
system. It should be independent of L. Figure 6 shows the average DOS at 
E = 0.5 £ for different system sizes. There are still some fluctuations present. 
These can in principle be reduced by using larger system sizes and increasing 
the number of samples. The fluctuations will be particularly inconvenient 
when trying to compute <r(E). 

The reduction of the DOS with increasing disorder strength can be un- 
derstood from a simple argument. If the DOS were constant for all energies 



12 



Alexander Croy, Rudolf A. Romer, and Michael Schreiber 




its value would be given by the inverse of the band width. In the Anderson 
model with box distribution for the on-site energies the band width increases 
linearly with the disorder strength W. The DOS in the Anderson model is not 
constant as a function of energy, nevertheless let us assume that for energies 
in the vicinity of the band centre the exact shape of the tails is not important. 
Therefore, 

^ - ^ , (29) 

shows a decrease of the DOS with W. Here B is an effective band width taking 
into account that the DOS is not a constant even for W—0. The parameter a 
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Fig. 6. Density of states vs disorder strength for E = 0.5 1 and L = 4, 6, 8, 10, 12, 14. 
Errors of one standard deviation are obtained from the ensemble average and shown 
for every 4th data point only. The solid line shows a fit to (29) for L = 6 to illustrate 
the reduction of p with increasing disorder strength. 

allows for deviations due to the shape of the tails. In Fig. 6, we show that the 
data are indeed well described by (29). 

6 Influence of the Metallic Leads 

As mentioned in the introduction most numerical studies of the conductance 
have been focused on the disorder transition at or in the vicinity of the band 
centre. Let us now set E = — 5t and calculate the conductance averages as 
before. The results for the typical conductance are shown in Fig. 7a. Ear- 
lier studies of the localization length provided evidence of a phase transition 
around W — 16.3 1 although the accuracy of the data was relatively poor [23]. 
Surprisingly, in Fig. 7a there seems to be no evidence of any transition nor of 
any systematic size dependence. The order of magnitude is also much smaller 
than in the case of E = 0.5 1, although one expects the conductance at the 
MIT to be roughly similar. 

The origin of this reduction can be understood from Fig. 8, which shows the 
DOS of a disordered sample and a clean system (i.e. without impurities and 
therefore without disorder) such as in the metallic leads. As already pointed 
out in Rcf. [19], the difference between the DOS in the leads and in the 
disordered region may lead to false results for the transport properties. Put 
to an extreme, if there are no states available at a certain energy in the 
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Fig. 7. System size dependence of the logarithm of the typical conductance for fixed 
energy E — —5t. Errors of one standard deviation are obtained from the ensemble 
average and are smaller than the symbol sizes. The lines are guides to the eye only. 
The upper plot was calculated using the metallic leads "as they are", i.e. the band 
centre of the leads coincides with the band centre in the disordered region. In the 
lower plot the band centre of the leads was "shifted" to the respective Fermi energy. 



leads, e.g. for \E\ > 6t, there will be no transport regardless of the DOS 
and the conductance in the disordered system at that energy. The DOS of 
the latter system becomes always broadened by the disorder. Therefore, using 
the standard setup of system and leads, it appears problematic to investigate 
transport properties at energies outside the ordered band. Additionally, for 
energies 3t ^ \E\ < 6t the DOS of the clean system is smaller than the 
disorder broadened DOS. Thus the transport properties that crucially depend 
on the DOS might be also changed in that energy range. The problems can 
be overcome by shifting the energy of the disordered region while keeping 
the Fermi energy in the leads in the lead-band centre (or vice versa). This 
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E/t 

Fig. 8. DOS of a clean system (full black line) and a disordered system (grey) with 
W = 12t and L = 21, obtained from diagonalising the Hamiltonian (1). The dashed 
lines indicate the band edges of the ordered system. 

is somewhat equivalent in spirit to applying a gate voltage to the disordered 
region and sweeping it — a technique similar to MOSFET experiments. The 
results for the typical conductance using this method are shown in Fig. 7. One 
can see some indication of scaling behaviour and also the order of magnitude 
is found to be comparable to the case of E = 0.5 1. Another possibility of 
avoiding the DOS mismatch is choosing a larger hopping parameter in the 
leads [45], which results in a larger bandwidth, but also a lower DOS. 

7 The MIT outside the band centre 

Knowing the difficulties involving the metallic leads and using the "shifting 
technique" explained in the last section, we now turn our attention to the less- 
studied problem of the MIT at fixed disorder. We set the disorder strength to 
W = 12t and again impose hard wall boundary conditions in the transverse 
direction [14]. We expect E c ss 8t from the earlier studies of the localization 
length [23]. Analogous to the transition for varying W we generate for each 
combination of Fermi energy and system size an ensemble of 10000 samples 
(except for L = 19 and L = 21, where 4000 and 2000 samples, respectively, 
were generated) and examine the energy and size dependence of the average 
and the typical conductance, (g^) and exp(ln<74), respectively. 

7.1 Energy Dependence of the DOS 

Before looking at the scaling behaviour of the conductance we have to make 
sure that the "shifting technique" indeed gives the right DOS outside the 
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ordered band. Additionally, we have to check the average DOS for being inde- 
pendent of the system size. In Fig. 9 we show the DOS obtained from diagonal- 
isation of the Anderson Hamiltonian with 30 configurations (using standard 
LAPACK subroutines) and the Green's function calculations. The Green's 
function data agree very well with the diagonalisation results, although there 
are still bumps around E = — 8.2 1 for small DOS values. The average DOS 
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Fig. 9. DOS vs energy for W = 12 1 obtained from the recursive Green's function 
method (circles with error bars obtained from the sample average) and from diag- 
onalising the Anderson Hamiltonian (histogram), for L A — 21 A = 9261. Also shown 
is a smoothed DOS (dashed line) obtained from the diagonalisation data using a 
Bezier spline. 

for different system sizes is shown in Fig. 10. For large energies the DOS is 
nearly independent of the system size. However, close to the band edge one 
can see fluctuations because in the tails there are only few states and thus 
many more samples are necessary to obtain a smooth DOS. 

7.2 Scaling Behaviour of the Conductance 

The size dependence of the average and the typical conductance is shown in 
Fig. 11. We find that for E/t < —8.2 the typical conductance is proportional 
to the system size L and the constant of proportionality is negative. This cor- 
responds to an exponential decay of the conductance with increasing L and is 
characteristic for insulating behaviour. Moreover, the constant of proportion- 
ality is the localisation length £. We find that £(E) diverges at some energy, 
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Fig. 10. DOS vs energy for different system sizes and W = 12 1 calculated with 
the recursive Green's function method. The data are averaged over 10000 disorder 
configurations except for L = 21 when 2000 samples have been used. The lines are 
guides to the eye only. Error bars are obtained from the sample average. 

which indicates a phase transition. This energy dependence of £ is shown in 
Fig. 12. 

For E/t > —8.05, (.94) is proportional to L. This indicates the metallic 
regime and the slope of (54} vs L is related to the d.c. conductivity. We 
fit the data in the respective regimes to the standard scaling form (26). The 
results for the critical exponent and the mobility edge are given in Table 2. 
The obtained values from both averages, (34) and (In 34), are consistent. The 
average value of v = 1.59 ± 0.18 is in agreement with results for conductance 
scaling at E/t — 0.5 and transfer-matrix calculations [9, 14]. 

average Emin/t £ max /t tir twr u E c /t 

JgZ) ^7A 3 2 1.60 ±0.18 -8.14 ±0.02 

(In g 4 ) -8.8 -7.85 3 2 1.58 ± 0.06 -8.185 ± 0.012 

Table 2. Best-fit estimates of the critical exponent and the mobility edge for both 
averages of qa using (26) with 71r = mi = 0. The system sizes used are L = 
11,13, 15,17, 19,21. 
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Fig. 11. System size dependence of the 4-point conductance averages (#4) and (In gn) 
for W = 12t and Fermi energies as given in the legend. Error bars are obtained from 
the ensemble average and shown for every second L. The dashed lines in the metallic 
regime indicate the fit result to (30) using the parameters of Table 2. In the insulating 
regime, linear functions for (In 54) = — L/£ + c have been used for fitting. 



7.3 Calculation of the D.C. Conductivity 

Let us now compute the d.c. conductivity from the conductance (g^E.L)). 
From Ohm's law, one naively expects the macroscopic conductivity to be the 
ratio of (g^E.L)) and L. There are, however, several complications. First, 
the mechanism of weak localisation gives rise to corrections to the classical 
behaviour for g 4 3> 1. Second, it is a priori not known if this relation still 
holds in the critical regime. And third, the expansion (27) does not yield a 
behaviour of the form 34 ae". 

In order to check our data for consistence with the anticipated power law 
(2) for the conductivity <j(E) in the critical regime, we assume the following 
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Fig. 12. Localisation length (for E < E c ) and correlation length (for E > _E C ) 
vs energy obtained from a linear fit to (In <? 4 ) = =F-k/£ + const., respectively. The 
error bars close to the transition have been truncated (arrows), because they extend 
beyond the plot boundaries. 



scaling law for the conductance, 

(34) = f{x v L) , (30) 

which results from setting b = x~ v m (24). Due to the relatively large error 
bars of (54} at the MIT as shown in Fig. 11, we might as well neglect the 
irrelevant scaling variable. Then we expand / as a Taylor series up to order 
71r and x m terms of e up to order mR in analogy to (27b) and (28). The best 
fit to our data is determined by minimising the \ 2 statistic. Using 7jr, = 3 
and wir = 2 we obtain for the critical values, v = 1.58 ± 0.18 and E c /t = 
—8.12 ± 0.03. These values are consistent with our previous fits. The linear 
term of the expansion of / corresponds to the conductivity close to the MIT. 
To estimate the quality of this procedure we also calculate the conductivity 
from the slope of a linear fit to (34) throughout the metallic regime, and from 
the ratio (54) / 'L as well. 

The resulting estimates of the conductivity are shown in Fig. 13. We find 
that the power law is in good agreement with the conductivity obtained from 
the linear fit to ((74) = oL + const, for E < —It. In this range it is also con- 
sistent with the ratio of (54) and L for the largest system computed (L = 21). 
Deviations occur for energies close to the MIT and for E > —It. In the criti- 
cal regime one can argue that for finite systems the conductance will always 
be larger than zero in the insulating regime because the localisation length 
becomes eventually larger than the system size. 
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Fig. 13. Conductivity a vs energy computed from (g±) / L for L — 21 (□), a linear 
fit with (34) = aL + const. (•) and a fit to the scaling law (30) (solid line). The 
dashed line indicates E c /t — —8.12. Error bars of (gA)/L represent the error-of- 
mean obtained from an ensemble average and are shown for every third E value. 
The dashed line indicates the position of E c and the inset shows the region close to 
E c in more detail. 



8 Conclusions 

We computed the conductance at T = and the DOS of the 3D Ander- 
son model of localisation. These properties were obtained from the recursive 
Green's function method in which semi-infinite metallic leads at both ends of 
the system were taken into account. 

We demonstrated how the difference in the DOS between the disordered 
region and the metallic leads has a significant influence on the results for the 
electronic properties at energies outside the band centre. This poses a big 
problem for the investigation of the MIT outside the band centre. We showed 
that by shifting the energy levels in the disordered region the mismatch can 
be reduced. In this case the average conductance and the typical conductance 
were found to be consistent with the one-parameter scaling theory at the 
transition at E c ^ 0. Using a finite-size-scaling analysis of the energy depen- 
dence of both conductance averages we obtained an average critical exponent 
v = 1.59 ± 0.18, which is in accordance with results for conductance scaling 
at E/t = 0.5, transfer-matrix calculations [9,11,14,26] and diagonalisation 
studies [10,12]. However, a thorough investigation of the influence of the leads 
is still lacking. It would also be interesting to see if these effects can be related 
to studies of ID multichannel systems with impurities [46] . 
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We calculated the d.c. conductivity from the system-size dependence of 
the average conductance and found it consistent with a power-law form at 
the MIT [47]. This strongly supports previous analytical and numerical cal- 
culations of thermoelectric properties reviewed in Ref. [3]. Similar results 
for topologically disordered Anderson models [48-50], random-hopping mod- 
els [37,51-53] and the interplay of disorder and many-body interaction [54-57] 
have been reported elsewhere. 
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